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Abstract 

We calculate the momentum dependent form factors of M — > 7*7* (M = vr°,ry) within the light- 
front quark model. Using the form factors, we examine the decays of M — > l + l~ , M — > 
and M — > l + l~l + l~ (l = e or u) and compare our results with the experimental data and other 
theoretical predictions. In particular, for tt° — > e + e~, we find that the decay branching ratio is 
6.68 x 10~ 8 , which is closed to the recent measurement of (7.48 ± 0.29 ± 0.25) x 10^ 8 by E799 of 
KTeV/Fermilab. 



I. INTRODUCTION 



The neutral pseudoscalar meson decays of M — > l + l~ , in particular Ki — > ft + fi~, have 
played very important roles to understand the Standard Model (SM). For the light pseu- 
doscalar mesons of 7r° and rj, the decays are dominated by the long distance (LD) contri- 
butions, described by the two photon intermediate state at the lowest order of QED. Since 
the short distance (SD) contributions in the SM are many orders of magnitude smaller, they 
can be neglected. Therefore, these decay modes are good processes to explore new physics 
beyend the SM. 

The measurement on this process by the KTeV-E799 experiment at Fermilab has given [lj 



B(n° e+e", x D > 0.95) = (6.44 ± 0.25 ± 0.22) x 10" 



(1) 



where xd = {m2 e /iTLn) 2 is the Dalitz variable with rri2e being the e + e~ mass. By extrapolat- 
ing the Dalitz branching ratio to the full range of xd with the overall radiative correction, 
one gets 



K°4e+e- = ( 7 ' 48 ± °- 29 ± °- 25 ) X 10 • 



(2) 



The decay of tt° — > e + e has been well studied theoretically over the years. However, the 



KTeV result in Eq. 



standard deviations |2|-|7| 



2]) disagrees with the some theoretical predictions about 1.5 ~ 3.3 



At the lowest order of QED, the decay branching ratio of tt° — > e + e is found to be[8MlO|: 
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where (3 = ^1 — ^m\jm\ and | A(mD\ 2 can be generally decomposed into |Im A{m^)\ 2 + 
I Re A{m^)\ 2 . Here, Im^4 denotes the absorptive contribution from the real photon in the 
intermediate state, which can be determined in a model-independent form 
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leading to the unitary bound on the branching ratio as 
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4.75 x 10" 



(5) 



The real part Re A is given by the dispersive one, which can be written as the sum of SD 
and LD contributions, 

Re A = Re A SD + Re A LD . (6) 

In the SM, the SD part is given by one-loop box and penguin diagrams Q]- The LD one 
involves the form factor related to the 7r°77 vertex. Using the form factor, the LD amplitude 
one has 

A J]_ [PV-(P-q) 2 ] F{q\{P- q f) 

* LD rr 2 m 2 J q q 2(P- q y[( q - Pe )2-ml} F(0,0) ' {) 

where P andp e are the pion and electron monenta, respectively. The function F(q 2 , (P—q) 2 ) 

is the double form factor of tt° — > 7*7*. This form factor contains the nontrivial dynamics 
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161 ] . In this paper, we 



of the process and has been studied in various models 0, 5^7, 
calculate the form factor F(q 2 , (P — q) 2 ) within the light-front quark model (LFQM) and 
use this form factor to evaluate the decays of tt° — » e + e _ and e + e~7. We will also study rj 
decays, which contain a dilepton or dileptons. 

This paper is organized as follows: In Sec. II, we present the relevant formulas for the 
matrix elements and form factors for M — > 7*7* (M = 7r°, 77). In Sec. Ill, we show our 
numerical results on the form factors and the branching ratios of meson M decays with 
dilepton. We give our conclusions in Sec. IV. 

II. THE FORM FACTORS 

To calculate M — > 7*7* (M = n°, rf) transition from factors within the LFQM, we have to 

xplicitly, 7T° may be described as (mm- dd)/V2 



decompose the mesons into QQ Fock states, 
and the valence state of rj can be written as jl7| 

\rj) = $ 8 cos8 P \uu + dd - 2ss)/VQ- $ 1 sin 9 P \uu + dd + ss)/V3 , (8) 

where $ 1,8 are the wave functions of the Fock states and 9p ~ —20° is the mixing angle. In 
the scheme of the QQ state, the amplitude of M — > 7*7* with CP conservation is given by: 

A(QQ(P) 7 *( gi , e a ) 1 *(q 2 , e 2 )) = ie 2 F QQ {q{, q 2 2 ) e^ pa 4 e\ q{ q^ , (9) 

where Fqq(q 2 ,q%) in Eq. d^J) is a symmetric function under the interchange of q\ and q\. 
From the quark- meson diagram depicted in Fig. 1, we get 




FIG. 1. Loop diagrams that contribute of ir c 
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+(Pi(3) -H-P3(i), Hmg), (10) 

where iV c is the number of colors and Ap is a vertex function which related to the QQ 
meson. In the light front (LF) approach, the LF meson wave function can be expressed by 
an anti-quark Q and a quark Q with the total momentum P as: 

\M(P,S,S,)) = J2 j [dpi][d P 2}2(2nf5 3 (P- Pl - p 2 ) 

Ai A2 

x$ s J'(z, k ± )b±(p 1 ,X 1 )d+(p 2 ,X 2 )\0), 

and 

dp + d 2 p± 



(11) 
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2(2vr) 3 ' 

where $^ A2 is the amplitude of the corresponding Q(Q) and p^ 2 ) is the on-mass shell LF 
momentum of the internal quark. In the momentum space, the wave function $^f z is given 
by 

^(k h k 2 ,X h X 2 ) = R s x %(z,k ± ) <P(z,k ± ), (13) 

where <f)(z, k±) describes the momentum distribution amplitude of the constituents in the 
bound state and R\A 2 constructs a spin state (S, S z ) out of light front helicity eigenstates 
(A1A2) 19]. The LF relative momentum variables (z, k±) are defined by 

z p+, p + = (1 - Z )P + , 



Pi 



Pix = zP± - kj 



P2± 



;i - z)p x + k A 



(14) 



The normalization condition of the meson state is given by 



(M(P', S', S' Z )\M(P, S, S z )) = 2(2vr) J P+8 s (P' - P)6 S , S 5 S , S , 



(15) 



which leads the momentum distribution amplitude <f)(z, k±) to 

dz d 2 k 



2(27ry 



(16) 



We note that Eq. ( TT3l) can, in fact, be expressed as a covariant form 2CH22] 

vtvi 
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z 1 — z 

In principle, the momentum distribution amplitude <f)(z, k±) can be obtained by solving the 
light-front QCD bound state equation 22] . However, before such first-principle solutions are 



available, we would have to be contented with phenomeno 



that has been used is the Gaussian type wave function 



ogical amplitudes. One example 



2J425J: 
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where N = 4(71-/0^)4, k = (k±, k z ), and k z defined through 
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M = E 1 + E 2 . 



(20) 



and dk z /dz = E x E 2 j z{\ — z)M . After integrating over p 3 in Eq. (fTU]) . we obtain 



A(QQ -> 7 *(?i) T*fe)) = e Q e Q N c p dp+ J 



d 2 p. 
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/ = Tr[ 75 (- ^3 + m ) ^Qfc + m Q ) + m Q )\ , V ~ ian = ^ + / L (22) 
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where the subscript {on} represents the on-shell partic 



function Ap from Eqs. ffTOj) . f fTTj) and f l2T|) . given by 
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es. One can extracted the vertex 



20 



2l|: 



AP - ^ W, (23) 



^ Plon P3on V2[M 2 - (mg-mg)' 

To calculated the trace /, we have used the definitions of the LF momentum variables 
(z(x), k±(k' ± )) and taken the frame with the transverse monentum (P — q 2 )± = for the QQ 
state(P) and photon(g 2 ) in Fig. la. Hence, the relevant quark variables are: 

P t=zP + , P t = (l-z)P + , Pl ^=zP L -k ± , P 3± = (l-z)P ± + k L . 

vt = xqf, pt = (1 - x)q£, p 2 ± = xq 2± - k' ± , p 31 _ = (1 - x)q 2± + k' ± . (24) 
At the quark loop, it requires that 

k± = (z - x)q 2± + k ± . (25) 

The trace I in Eq. (|22|) can be easily carried out. Thus, the form factor F(qf, q 2 ) in Eq. (Q 
can be found to be: 

where cq is the quark electric charge factor and 



2(1 — z) / d/c, / 
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If gi and g 2 are on mass shell where r = 1, the form factors of tt — > 77 and 77—7-77 can be 
written as 

iV c /" dxd 2 k± $ (x, fcj_) f4 m u 1 m<f 



F^ 77 (0,0) = 8v / 2y^ y 



2(2vr) 3 1-z \9m£ + fci 9m 2 + / ,2 j ' 

p rnm-lfi /S f dx( ^ k ± j ® 8 {x,k 2 ± )cos6 P (4 m u 1 m d 2 m s \ 
W ' J V 3 7 2(2tt) 3 I (l-x)V6 V9m2 + fci + 9m2 + fc2 9 + ^ y< 
$ : (x, /c 2 ) sin 6*p / 4 m M 1 1 m s 



l-x)V3 \9m 2 u + kl 9m 2 + kl 9m 2 + k] 



28 
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III. NUMERICAL RESULT 



To numerically calculate the transition form factors of 7r° and 77 in Eq. (l26|) and ( )28|) . we 
need to specify the parameters appearing in <p(x,k±). To constrain the quark masses of 
m UjdjS and the meson scale parameters of um in E q. (I2"6"j) . we use the meson decay constants 
fu and its branching ratios of M — > 27, given Toy [271] 

/>=132MeV, /®=169MeV, /,}=145MeV. (29) 

and 



J Br 7r o_, 27 



(98.832 ±0.034)%, Br 



(39.30 ±0.2) 



respectively. Here, the explicit expression of ju is given by 

. \fW r f dx d 2 k 



26] 
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m d P \F(0,0) 



(30) 



(31) 



(32) 



64vrr P 

we find that |-F(0, 0) 7r o( r) )^ 2 7l = 0.274(0.272) in GeV^ 1 . As an illustration, we extracte 
m u = rrid = 0.24, m s = 0.38 and u n = 0.33, u v i = 0.42, u v $ = 0.58 in GeV, which will be 
used in our following numerical calculations. 



e T e 7 



We now examine process of tt° — > e + e~ / y with the form factor in EqJ26]). The interaction 
between the photon and leptons is given by the conventional QED 14j,|28j. One easily obtains 
the differential decay rate 



dT(n° e + e~j) _ 2a 1 
r(vr° 77) dq\ ~ 3^qf ' ~ 



1/2 



, 2 m 
1 + —T- 



\f(t)\\ (33) 



where fit) = F-^^q 2 ,0)/ F n (0,0) and t = q\jm\. Obviously, the branching ratio of n° — > 
e + e~7 in the Eq.( l33|) depends on the factor of l/qf- The function of f(t) is an analytic 
function in the entire physics region of Am 2 < q\ < m 2 , related to 



± 



x(l — x)qf — m\ — k\ m\ + k 



- -(m u o m d ) 



(34) 
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Integrating over q\ in Eq. ( 133|) . we obtain the branching ratio 

r ff^*7> = 1.18X10-', (35) 
r(7r u — > 77 

which agrees well with those by QEd[i^,[i5] and vector meson dominance (VMD) model[l^]. 
Our result is also close the experimental data: <B^o^ e+e - 7 = (1.198 ± 0.032) x 10~ 2 [27| . 

Similarly, the branching ratios of 77 — > e + e~7 and 77 — > /i + /i~7 which normalized with 17 
tatal width are found to be 

T(i] — > e + e~7) 
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r(r/^/i+/i 7) 



6.95 x 10- 



2.94 x 10" 



(36) 



Ours result of 77 — > e + e 7 is smaller than that in the CLEO data 291 but larger than the one 

n n 

in Ref.[31J. However, for the mode of 77 -± u + /i 7, our result agrees with Ref.[31J as well as 



that by the effective mass theory (EMT) 321 ] . Furthermore, our predictions in the two decay 

33] and the PDgQ. 



modes agree well with the experimental data in CELSIUS 



B. 7T° e+e-e+e" and 77 -> 1+ 1~ i+ 1~ (£ = e, 77) 

We examine the double lepton-pair decay of 7r° — » e + e~e + e~ with the form factors in 
Eq. (I2"B"1) . The decay matrix element is calculated by the conventional QED with the inter- 
action of 7T° and two photons and the differential decay rate is given by 



dr(vr -> e+e-e+e" 



r(7T° 77) dg 2 dq\ q\q\ V3tt 



a 



^(0,0) 



A 3/2 fl,4'4) GKg?)C,(? 2 2 ).(37) 



where 



A(a, 6, c) = a 2 + 6 2 + c 2 — 2(ab + be + ca) 



Gi{q 2 



4 mi 



2\ 1/2 
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2 mi 



After the integrations over qf and gf, we obtain the branching ratio as follows: 



7r u — >e T e _ e" r e" 



r(vr° -)■ e+e-e+e- 

r(7r° — >■ 77) 



3.29 x 10" 



(38) 



(39) 



which is smaller than that in Ref. 14j . but larger than the one in Ref. 15|] slightly. However, 
all results are consistent with the experimental data. We note that even if the form factor 
is replaced by an on-shell constant with F(g 2 ,g 2 ) = F(0,0), the branching ratio is found 
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to be very close to the result in Eq. ( 139]) . It might be a good approximation to neglect the 
momentum dependence of the form factor for the decay. 

We can also perform the similar calculations for r\ — > l + l~l + l~{l = e or ji) and we find 



B^ e + e - e + e - = 2.47 X 10 



B^e+e-At+M" = 5 ' 83 X 10 



B^+m-^+m- = 1.68 x 10- 9 . (40) 

Our result on B v ^. e + e - e + e - is in good agreement with the experimental data £T^ e+e _ e+e _ = 
(2.7^2^ stat ± O.lsyst) x 10~ 5 33| and Ref . [31I . For other modes, currently, our theoretical 
predictions are many orders of magnitude smaller than the experimental upper bounds 
33j. 



27|, 



c. 7r°(r]) -> e+r 

We first calculate the real part of Re Ald in Eq- © at the pion momentum limit of 
P 2 0. At this limit, the relevant form factor of Eq. ( 126]) . given by a triangular quark loop, 
would be simplify to 



2 (2tt)" 1 -' 1 - z 

m„ 1 m d 



9 sc(l — x)g 2 — m 2 — fc 2 9 x(l — x)g 2 — m 2 d — k 



(41) 



One could easily find 

Re A LD (0) ~ -20.74 . (42) 

The numerical result is in agreement with the most vector meson dominance (VMD) model 
at P 2 — > 0. This implies the equivalence between the VMD and LFQM descriptions on the 
form factors of hadrons with the relevant vector meson mass of My ~ 2m M in the VMD. 
To illustrate Re Ald{q 2 ) m the range — m 2 > q 2 > m 2 , we use the dispersive framework 



30] 



proposed in Ref. j^]. The real part may be written by a once-subtracted dispersion relation 



a 2 r°° Im A(a' 2 ) 
Re A LD {q 2 ) = Re ,4(0) + ±- \ dq' 2 a (43) 

Jo {q' 2 -q 2 )q' 2 



Extrapolating from q 2 = to m%, we find Re Ald{j^) = 11.18. Since the SD part of Re Asd 
can be neglected, we get the branching ratio of the real part in Eq.(l) to be 1.93 x 10~ 8 . 
The total decay branching ratio is about 6.68 x 10 -8 . Our prediction is smaller than the 
experimental value of BjS^e+e- = ( 7 - 48 ± °' 29 ± °' 25 ) x 10_8 measured b Y KTeV. We note 
that our result is larger than the values of (6.41 ± 0.19) x 10~ 8 and 6 x 10~ 8 calculated in 
Ref. with the VMD and quark model(QM), respectively, but closed to (7± 1) x 10 8 in 
the Chiral Perturbation Theory(ChPT)|a|. It is clear that we provide a method to calculate 
the form factor of tt° — > 7*7* and get a result in 7r° — > e + e~ within the LFQM. 

The i] — > decay can be analyzed in a similar technique as 7r° — > e + e~ . In the 

momentum limit P 2 — > 0, we obtained 

Re A { 2e)L D (0) ^ -22 A3 , 

Re A<M LD (0) ^ -6.48 . (44) 
Form the dispersive integral in Eq. (j4"3]) and Eq. fj44l) . one obtains 

Re A { 2e)L D (m*) ~ 27.11 , 

Re A {2fl)LD (m 2 ) ~ -2.81 . (45) 

The SD contributions to the decays can be still ignored and the total branching ratios are 
given by 

B n ^ e+e - = 4.47 x 10- 9 , 

B^m+m- = 5-47 x 10- 6 . (46) 

One notes that the value of B v ~> e + e - is larger than the CLEO result {2}. For the mode of 
rj — > it is consistent with the CLEO [2] and VMD results (3]. It also agrees with the 

PDG data of 5.8 ± 0.8 x 10~ 5 . 

We summarized the related experimental and theoretical values of the decay branching 
ratios of 7r° — > e + e~7, 7r° — > e + e~e + e~ and 7r° — > e + e~ in Table I and r] — > / + / _ 7, r] — > 
l+l~l + l- and 77 ->• l+r in Table II. 



IV. CONCLUSIONS 
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TABLE I. Summary of the decays of tt° with lepton pair. 



Br 


Exp. data 


This work 


Other models 


10 2 £> e + e - 7 


1.174 ±0.035 [22] 


1.18 


1.18[14J[15][16] 


10 5 B e + e - e + e - 


3.34 ± 0.16[27] 


3.29 


3.28[14], 3.46[15] 


10 8 B e+e - 


7.48 ± 0.29 ± 0.25[1, 2] 
6.46 ± 0.33[27] 


6.68 


7±1[3], 8.3±0.4[4], 6.41 ±0.19[5j, 6[6J, 
< 4.7[7], 6.23 ± 0.09[2, 30J 



TABLE II. Summary of the decays of rj with lepton pair. 



Br 


Exp. data 


This work 


Other models 


io 3 iw 7 


7.8±0.5 sia i ±0.7 syst [33j 
7.0±0.7[27] 


6.95 


9.4±0.7[29] , 
6.31 - 6.46[31], 6.5[32] 


10 4 B M +M-7 


3.1±0.4[27J 


6.95 


2.14 - 3.01[31], 3.0[32] 


10 5 £> e + e + e - e - 


2 - 7 i|7^at =t 0.1 s?/s t[33j 
< 6.9[27J 


2.47 


2.49 - 2.62[31] 


!0 7 S M + M - e + e - 


< 1.6 x 10 3 [27] 


5.83 


1.57 -2.21 [31] 




< 3.6 x 10 5 [27J 


1.68 






< 2.7 x 10 4 [27] 


4.47 


13.7[5] , 4.60 ±0.06 [2, 30] 


10 6 Vr 


5.8±0.8[27] 


5.47 


5.8±0.2[3], 11.4[5] 
5.11 ±0.20[2, 30], 5.2 ± 1.2[34] 



We have calculated the form factors of P ->. 7*7* (P = vr°, T]) directly within the LFQM. In 
our calculations, we have adopted the Gaussian-type wave function and evaluated the form 
factors for the momentum dependences in the energy regions from q 2 = to m 2 P . Using the 
form factors, we have examined 7r° — > e + e _ 7 and 7r° — > e + e~e + e~ and shown that our results 
on the decay branching ratios agree well with the experimental data shown in Table. I. Our 



predicted values are also close to those in the QED and VMD models|14H10j|. For 7r u — > e^e" 
we have found that S„-o^. e + e - is 6.68 x 10~ 8 , which agrees with (7±1) x 10~ 8 in the ChPT 
but larger than those in Refs.p-Li]]. We have demonstrated that the long-distance dispersive 
contribution in this model is possibly small. However, like other theoretical predictions, our 
result for 7r° — > e + e~ is also slightly smaller than the experimental data. Clearly, further 
theoretical studies as well as more precise experimental data such as those from the KTeV- 
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E799 experiment at Fermilab on the spectra of the decays with lepton pair are needed. 
About the rj decays, our results are all consistent with the experimental data. In particular, 
the branching ratios of rj — > 2e2/i, rj — > Afi and rj — > 2e are expected to be 4~5 orders of 
magnitude lower than the current experimental upper limits. 
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